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ABSTRACT 


The CH-46 Tie Bar is a multiple leaf, stainless steel system that attaches the rotor blade to 
the rotating hub and carries rotor blade centrifugal force. The Tie Bar twists as the rotor blade 
feathers and bends slightly as the Pitch Housing bends. Onginally designed to last 3,000 hours 
without any component breakages, the Tie Bar has experienced field failures that have necessitated 
inspections every LO flight hours. Traditional fatigue qualification and life methodology have not 
reconciled the unusually large number of CH-46 Tie Bar component failures. One CH-46 
experienced 15, 10, and 6 component failures on the three aft rotor head Tie Bars after only 93 


flight hours. One complete Tie Bar (140 components) failed and caused an aircraft crash. 


This research developed a probability modeling predictive method that can aid the design 
of a Tie Bar that will meet desired life and reliability levels. Laboratory coupon tests were 
interpreted to provide material properties for several manufacturing processes. Finite element 
analysis of the current Tie Bar design as well as two modified designs was used to determine the 
Tie Bar stress state for many component failure combinations. The reduced constraint design 
produced a dramatic stress concentration reduction and may provide large gains in life over the 


current and elliptical slot designs. 
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I. INTRODUCTION 


A. CH-46 HELICOPTER 

The CH-46 helicopter, pictured below in Figure I-1 and Figure 1-2, is operated by both 
the U.S. Navy and the U.S, Marine Corps. Each rotary wing head on the CH-46 has three rotor 
blades. A diagram of the rotary wing head is shown in Figure [-3. This investigation addresses the 


analysis and prediction of field failures of the Tie Bars and options for Tie Bar redesign. 


GROUND 
LINE 


Figure I-1 CH-46 Helicopter, Side View 


Figure I-2 CH-46 Helicopter, Top View 
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Figure I-3 CH-46 Rotary Wing Head 


B. PITCH ASSEMBLY 


The exploded view of the CH-46 Rotary Wing Head pitch assembly, Figure I-4, depicts 
the Tie Bar and the other rotor blade pitch change components. The cutaway view, Figure I-5, 
illustrates the Tie Bar’s location inside the Pitch Housing and Pitch Shaft. 

The Tie Bar carries the rotor blade centrifugal force. One end of the Tie Bar is connected 
to the Pitch Housing; the other end is connected to the Pitch Shaft. There are roller bearings 
between the pitch shaft and pitch housing that allow relative rotation and sliding between the two 
parts. Since the Tie Bar is connected to both the Pitch Housing and the Pitch Shaft it twists as 
control forces are applied to the Pitch Housing to feather the rotor blade. To reduce control loads, 
the Tie Bar is designed to be soft in torsion. A flatwise bending moment due to flapping applies a 
bending moment to the pitch shaft and pitch housing which in tum bends the Tie Bar. 
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Figure I-4 CH-46 Rotary Wing Head Pitch Change Assembly, Exploded View 
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Figure I-5 CH-46 Rotary Wing Head Pitch Change Assembly, Cut-Away View 


C. TIE BAR 


The Tie Bar is constructed of 301 Stainless Steel. It consists of 35 slotted straps that are 
stacked with 301 Stainless Steel spacers in between. The slots, in essence, make each individual 
Tie Bar strap a redundant four component assembly. The “components” are lettered A through D 
in Figure I-6, The Tie Bar stack is held together by six hollow pins that are expanded by driving a 


metal ball through them dunng assembly. 
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Figure I-6 Tie Bar Strap Geometry 


Designed to last 3,000 hours without any component breakages, the Tie Bar is now 
inspected every 10 flight hours. Normal fatigue qualification and life methodology have been 
unsuccessful in predicting an unusually large number of CH-46 Tie Bar component failures. One 
CH-46 experienced 15, 10, and 6 component failures on the three aft rotor head Tie Bars after only 
93 flight hours. A complete Tie Bar failure has caused one aircraft crash. The cause for early in- 
service failure of Tie Bar elements has not been determined. Improved manufacturing methods and 


quality control have not eliminated the problem. 


D. PROBABILITY MODELING 


To predict the life of a part by statistics with a high confidence level for safety related 
components, the number of test specimens must be at least an order of magnitude higher than the 
desired reliability. For example, to obtain a reliability level of one failed part in 100,000 (10°), at 
least 1 million (10°*') parts need to be tested. Clearly, this is impractical. In addition, the statistical 
approach requires building the actual part before the reliability level can be determined. To solve 
the Tie Bar failure problem, a modeling method needs to be developed that can aid the design of a 


Tie Bar that will meet desired life and reliability levels. 


Since a statistical approach is unrealistic, a probability approach has been chosen. 
Through a combination of finite element modeling, coupon test analysis, fatigue modeling, and 


probability modeling, this investigation attempts to develop a model to predict Tie Bar life. 


1. COMBINATORIAL FAILURE SEQUENCES 


O'Connor [Ref. |] studied the Tie Bar problem in a previous investigation. He identified 
that the probability of failure of a system with multiple components could be determined in terms 
of the probability of failure of the components. That is, the probability of failure of the Tie Bar can 
be determined in terms of the probability of failure of the straps which can in tum be determined in 
terms of the probability of failure sequences of the components within the straps. He also identified 
that each component of a Tie Bar strap can be thought of as a chain in which failure of the weakest 


link results in failure of that entire component. This lead to adoption of a Weibull distribution for 


F@&) = ie -(5) . 


Equation II-1 


the individual strap components. 


Applied to the Tie Bar, 
Fle, (2)|=1- os Ga . 


Equation II-2 


where F’,, [x F (L)| is the probability of failure of component A within a strap under stress x, in 


the i” configuration (one, two, three, or four component failures), xa(L) is the applied stress on 
component A caused by the global load L, B is a scale parameter of the material strength, and « is 
a shape parameter of the maternal variability. It is important to note that the components have 
identical failure probability parameters (a, §) because the material is identical, but because the 
stress states are not the same, each component within a strap has a different probability of failure. 
Finite element analysis is used here to calculate the stress states for input into the probability of 
failure, Equation II-2. Because of the local stress dependence on the combinatorial failure 


sequences, the probability of failure of the Tie Bar is no longer Weibull. 


A total of 16 failure states exist for a four component Tie Bar strap. Failure states are 
numbered using a binary code. A I represents an unbroken Tie Bar strap component and a 0 
represents a broken component. Table II-1 shows the possible failure states and a decimal 


equivalent name for each failure state. 


Component State Decimal 
A B Cc Equivalent 
1S 
14 
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1 represents an unbroken Tie Bar strap component 


0 represents a broken component 


Table II-1 Tie Bar Strap Failure States 


O’Connor [Ref. |] developed an algorithm to calculate the probability of failure of a Tie 
Bar leaf given the initial configuration; zero, one, two, three, or four components failed and the 
final configuration one, two, three, four, or all failed. He first developed an algonthm to calculate a 
three component system and found that the local stress can be calculated for each break 
configuration, but the relative strengths cannot be determined unless the failure path which lead to 


that configuration is known. His three component “failure tree” is shown in Figure II-1. 


A=B=C Equally Weak 


B=C Weakest 


A=C Weakest 


A=B Weakest 


C Weakest 


XM Indicates 
Failure 
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A Weakest 


Figure II-1 Three Component Failure Tree 
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The algorithm to generate the probability of system failure given the component failure 
first requires calculating the probability of failure of a sequence of consecutive states, (one branch 
or path of the failure tree). The probability of failure of one sequence is the intersection or product 
of the transitional probabilities. The probability of failure of the system is the union or sum of the 
probability of failure of each of the sequences. The appendix of O’Connor’s work (Ref. 1] lists all 
the possible failure permutations for a four component system given that none or any one, two, or 


three of the four components has already failed. 


tf. FINITE ELEMENT MODELS 


In order to calculate the probability of failure of components of the strap, the stress state at 
the prospective failure location must be known for each of the physical states, (none failed and 
partial failures in Figure II-1). Finite element analysis is used to calculate such states for three 


different designs. These designs include: 
A. The current design 
B. Anelliptical slot design proposed by Boeing Vertol 


C. A reduced constraint design proposed in this investigation 


A. CURRENT DESIGN 

A Finite Element model of the CH-46 Tie Bar, Figure III-1, was created using 
MSC/NASTRAN for Windows. The model consists of the top and bottom straps of the Tie Bar, 
the three expanded pins at each end of the Tie Bar that hold the assembly together, and the two 
large pitch housing/pitch shaft pins (parts 162 and 167 of Figure I-4) that connect the Tie Bar to 
the pitch shaft and pitch housing. The straps are modeled using quad PLATE elements. The top 


and bottom straps are meshed identically. 


Current Tie-Bar 


Figure III-1 Current Tie Bar Finite Element Model 


The expanded assembly pins and pitch housing/pitch shaft pins are modeled using RIGID 
elements. The assembly pin RIGID elements tie all degrees of freedom of a bottom strap node to a 
corresponding top strap node. The pitch housing/pitch shaft pins (Figure III-2) are modeled by 
two vertical RIGID elements to which “bicycle spoke” RIGID elements are radially attached. 
Boundary conditions are applied to the center node of the pitch housing/pitch shaft pins to model 
the displacements of the Tie Bar in flight. 
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Figure III-2 Pitch housing/pitch shaft Pin Detail 


The pitch housing/pitch shaft pins in this finite element model are attached to the tie bar 
straps. These pins in the actual Tie Bar are not attached. Attaching the pins to the Tie Bar straps 
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provides expedient modeling of the equivalent boundary conditions. The local stress around the pin 
holes is affected by this expediency; however, service history has shown that failures occur in the 4 
components of each Tie Bar strap, not the areas around the holes. Saint-Venant’s principle assures 
that the stresses at actual failure locations, which are far away from the pin holes, are accurately 

modeled. By applying the boundary conditions to the attached model pins, all the displacement is 

applied directly to the 4 strap components. The pitch housing/pitch shaft pin holes remain circular; 
they do not elongate, Therefore, the four strap components experience slightly more strain than the 


real Tie Bar, and the model is conservative. 


The Tie Bar straps are positioned to simulate the maximum allowable Tie Bar stack 
thickness. The PLATE elements around the pitch housing/pitch shaft pins are thicker than the rest 
of the Tie Bar, It is assumed that the spacers on both sides of the outermost Tie Bar straps 


contribute to the strength and stiffness of that region. Figure III-3 illustrates the thickened region. 


Current Tie-Bar 


Figure II]-3 Side View of Tie Bar Model Showing Thickened Regions 


B. ELLIPTICAL SLOT DESIGN 


The Tie Bar manufacturer, Boeing Vertol, has suggested modifying the geometry of the 
Tie Bar strap slots to reduce stress concentration and increase Tie Bar life. The focus of the 
modification is to change the previously semi-circular slot end to a semi-elliptical slot end as shown 
in Figure III-4. This geometric modification was made to the finite element model to replicate the 


new geometry. 
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Figure III-4 Tie Bar Slot Geometry Modification 


C. REDUCED CONSTRAINT DESIGN 


In studying the kinematics of the Tie Bar, it was observed that rotor blade feathering 
causes straps away from the Tie Bar neutral axis to experience not only a twisting deformation but 
also a parallel shear. If the three swedged assembly pins around the pitch shaft/pitch housing pins 
are removed, the reduced constraint significantly minimizes the shear deformation. Without 
changing the geometry from the current Tie Bar design, another set of models was generated by 


removing the three swedged assembly pins from each end of the Tie Bar model. 
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Figure III-5 | Reduced Constraint Tie Bar Finite Element Model 


The straps for the reduced constraint models are only held together by the pitch 
housing/pitch shaft pin on either end. Ina real Tie Bar of this design, the straps could be held 
together by a sleeve with rolled-over ends, Figure III-6. The assembly pin holes would not be 


required, and the Tie Bar strap ends would be rounded off to maintain clearance inside the pitch 


Side view Top view 


Figure [II-6 Sleeve to Assemble Stack of Reduced Constraint Tie Bar 


shaft when the straps rotate. Also, a low friction material such as Teflon could be used for the 
spacers between straps to reduce friction. These models were easy to create from the current Tie 


Bar models and permitted direct comparison. 


Whereas the current and elliptical slot design models have all degrees of freedom (DOF) of 
the pitch housing/pitch shaft pin end nodes tied to the pitch housing/pitch shaft pin master node, 


1S 


the reduced constraint model allows the pitch housing/pitch shaft pin end nodes to rotate about the 
pin z-axis relative to the pitch housing/pitch shaft pin master node. This design, Figure III-7, 


alleviates in-plane shear and bending of the Tie Bar straps. 
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Figure III-7 Pitch housing/pitch shaft Pin with Rotating Ends 


D. FAILURE STATES 


Models with break sites were created from the unbroken, “parent” models. Individual top 
strap components were “broken” to simulate breakage in service. This was accomplished by 
disconnecting the nodes at the prospective failure sites to model the crack trajectory. With the 
nodes disconnected, the affected PLATE elements in the model were shrunk to show the “break” in 
the Tie Bar strap. 


As stated previously, a total of 16 failure states exists for a four component Tie Bar strap. 
The possible failure states are shown in Table II-1. Six of these states are eliminated using 
symmetry, and the all-broken state does not require analysis. So, eight models were “grown” from 
the original model to recreate the possible failure states. Figure III-8 illustrates a broken Tie Bar 


model. 


Figure III-8 Broken Tie Bar Strap 


E. BOUNDARY CONDITIONS 

The boundary conditions for the Tie Bar were supplied by the Tie Bar manufacturer as a 
flight spectrum. For each flight condition, values of centrifugal force, twist, and bending were 
provided. The flight conditions consist of G-A-G, Pre-Flight, HSLF, Maneuver, and Autorotation, 
and they are listed with their corresponding Tie Bar displacements in Table III- 1. 


0>113%0 
0298 RPM0 


Pre-Flight 
Maintenance 


HSLF 


Maneuver 


0>118%-0 
0312 RPM>0 


100% Nr 
264 RPM 
105% Nr 
277 RPM 
105% Nr 
277 RPM 
105% Nr 
277 RPM 


Limit 


Ultimate 


105% Nr 
277 RPM 
113% Nr 
298 RPM 
113% Nr 
298 RPM 
113% Nr 
298 RPM 
113% Nr 
298 RPM 
118% Nr 
312 RPM 
133% Nr 
351 RPM 


076250 Ib>0 
083150 Ib>0 
0 
0 
59715 Ib 
65835 Ib 
65835 Ib 
65835 Ib 


65835 lb 


76250 Ib 


76250 lb 


76250 lb 
83150 Ib 


105630 Ib 


139960 Ib 


-25°—- 16° +22° 
-25°- 16° +22° 
-25°-+22° 


+50° 
0+6° 


~11.6°+4.5° 


-10°+6° 


-6°+9.4° 


-11.6°+4.5° 


-20° 


-20° 


0>+0.365°0 
0+0.365°0 
Negligible 


0 
+0.23° 


* Twist cycle is once per rotor revolution for HSLF, Maneuvers, and Autorotation conditions. 
@ Relative angle of attachment pins referenced to the vertical pin. Bending cycle is two times per rotor 


revolution. 


Table III-1 


1. Flight Condition Description 


Tie Bar Loading Conditions 


6000 
5 


3.33x10’ 


5.34x10° 


0.55x 10° 


6000 


3000 
16700 
Once per 


Lifetime 
N/A 


The G-A-G condition is a ground-air-ground condition. The first half of the G-A-G cycle 


includes engaging the aircraft rotors, checking the rotor RPM at maximum throttle with the 


collective bottomed, and lifting the aircraft off the ground. The second half of the G-A-G cycle 


includes the aircraft landing and rotor shutdown. 
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Preflight checks occur on the ground with the rotors disengaged. The collective is moved 
from full down to full up. The cyclic is moved full-forward, full aft, full left, and full right during 


these checks. These control movements twist the Tie Bar. 


HSLF is honzontal steady level flight. It is a cruise condition without climbs, descents, or 


tums. 


Autorotation is a power off maneuver in which the collective is lowered and the aircraft's 
descent is slowed by the windmilling of the rotor blades. As shown in Table III-2, autorotation is 


generally flown at a higher rotor RPM than the other flight conditions. 


2. Rotor System Behavior 


O'Connor [Ref. |] conducted a kinematic analysis of a single rotor system in previous Tie 
Bar research to establish the phasing of the loads and the kinematics of the Tie Bar. A fully 
articulated system, the CH-46 rotor system has both lead-lag and flapping hinges to provide 
independent in-plane and out of plane motion of each rotor blade. The rotor blade flapping hinge is 
designed into the rotor system to prevent development of an undesired rolling moment on the 
helicopter. Flapping occurs because the advancing rotor blade sees a higher relative velocity and 
hence higher lift than the retreating blade. Maximum upward flapping occurs at y = 180°, and 
minimum flapping occurs at y = 0°. Through his research, O'Connor also discovered that the 
maximum flatwise bending moment occurs at y = 90°->120° (upward) and y = 270° (downward). 


The rotor angle convention is shown in Figure III-9. 


y = 90° y = 270° 


= 180° 


Flight Direction 


Figure III-9 Rotor Angle Convention 


Lead-lag is caused by cyclically varying drag and inertia loads. As a blade flaps up, its 
center of gravity moves inward. Since angular momentum must be conserved, the blade accelerates 
in plane to compensate for the decrease in moment of inertia. The change in inertia of this upward 
flapping blade is resisted by the other blades and an in-plane moment is produced. The lead-lag 
hinge allows the blade to accelerate forward. As the blade flaps down, the center of gravity moves 


out, and the blade decelerates. 


Maximum twist or feathering of the rotor blade occurs at y = 270° where the rotor blade 
increases pitch to compensate for the lower relative velocity. Minimum feathering occurs at y = 


90° where relative velocity is highest. 
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O'Connor [Ref. |] also researched the phasing of the Tie Bar bending and twist. He noted 
that the lead-lag produced negligible bending of the Tie Bar. His study of the rotor blade dynamics 


produced the results in Table III-2. 


Rotor Position - y Centrifugal Force Feathenng Bending 


(degrees) 
constant 
constant 


constant 


270 constant 


Table III-2 Rotor Phasing 


3. Conversion of Flight Conditions to Model Boundary Conditions 


Three primary load cases were run for each Tie Bar design and break state. The three load 


cases were: 
1. 0.2 inch extension at each pitch housing/pitch shaft pin. (X-translation); Figure III-10 
2. 0.2 radian twist at each pitch housing/pitch shaft pin. (X-rotation); Figure III-11 


3. 0.2 radian bending at each pitch housing/pitch shaft pin. (Y-rotation); Figure I-12 
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Outour Set: Extanmon 10.2 
DalormesiO. 1 Total Tanwevon 


Figure III-10 Primary Extension Load Case (Side View) 
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Figure III-11 Primary Twist Load Case 
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Figure IIf-12 Primary Bending Load Case 


These three basic load cases were scaled and linearly combined to recreate each flight 
condition at y = 0°, 90°, 180°, and 270°. The twist and bending boundary conditions were scaled 
directly. The extension boundary conditions were given as centrifugal forces instead of 
displacements. To scale these boundary conditions, the force to extend each unbroken model 0.2 
inches was read from the NASTRAN output. This value was divided by 2 to find the force 
required to deform | unbroken strap 0.2 inches. This result was multiplied by 35 to determine the 


force required to deform a whole 35-strap, unbroken Tie Bar 0.2 inches. 
(Fore = (Far), +2x35, 


Equation III-1 


where (F,2) w» is the force required to deform a whole Tie Bar 0.2 inches, and (F, al is the force 
required to deform a two-strap NASTRAN model. 


To calculate the extension boundary condition scaling factor the following equation was 


used: 
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CF 
(, 2)as - (Fy 2 dis) +34 *®, 2 Mis 


Extension_scale_ factor = 


Equation III-2 


where CF is the centrifugal force for the given flight condition, and (F, 7 I. is the force required to 


deform one unbroken strap 0.2 inches. 
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IV. FINITE ELEMENT MODEL RESULTS 


Service failures of the current Tie-Bar and finite element model stress contour plots 
indicate that the areas around the Tie-Bar slot ends are most cnitical. At every spatial location, the 
3 scalar components of the stress tensor can be combined to form a scalar value based on the 
failure process. Since stainless steel is a ductile material, the von Mises failure stress is used. The 
NASTRAN contour plots that follow, Figure I[V-1, Figure [V-2, and Figure [V-3 show the 
deformation and von Mises stress contour of the top strap of the three Tie-Bar designs when the 
Maneuver 4 flight condition is applied at y = 270. While the Maneuver 4 flight condition at y = 
270 loads the top Tie Bar strap more heavily than the bottom one, it qualitatively represents typical 
behavior for the models. The Current Tie Bar and the Elliptical Slot Tie Bar models exhibit in- 
plane bending caused by the assembly pins. The assembly pins constrain the rotation of the Tie Bar 
straps about the pitch shaft/pitch housing pins. The Reduced Constraint Tie Bar model exhibits no 


in-plane bending since the assembly pins have been removed. 


The current Tie-Bar model shows areas of high stress where the slots transition from a 
straight to a radius. The Elliptical Slot Tie-Bar model shows stress concentration in the same 
areas, but at a slightly lower value. The Reduced Constraint Tie-Bar Model does not show any of 


the stress concentration of the first two models around the slot ends. 


To begin understanding the effect of the new designs on the stress in the Tie-Bar 
quantitatively, several elements of interest were chosen. These elements are located at the transition 


between the straight part of the slots and the radiused ends. The element numbers are shown in the 


“element map”, Figure IV-4. 


Y 


Z—_ -x 


Output Set: Maneuver 4, Psi = 270 
Deformed(0, 131): Total Transtaticn 
Contour: Plate Top VonMises Stress 


$7809. 51750, 46000, 40Z50. MSH. 28750. 23000. 


(Note shear deformation, left side displaced toward Y+ (up) and right side displaced toward Y- (down)). 


Figure IV-1 


Current Tie-Bar, Top Strap, Maneuver 4, ‘¥=270 
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$3250. 57500, 51750. 46000. 


zZ—__-x 


Output Set: Maneuver 4, Psi=270 
Deformec{0.168): Total Transietion 
Contour: Plate Top VonMises Stress 


40286. 34500 


Figure 1V-3 


Reduced Constraint Tie- 


Bar, 


Top Strap, 


Maneuver 4, 


Y=270 
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Focusing on these “critical elements”, the effect of the design changes was first evaluated 
by looking at the “critical elements” with the primary load cases applied individually. This allowed 


study of individual effects on the Tie Bar behavior. 


For the extension only load case, the Tie-Bar was elongated 0.02 inches, a typical flight 
value. When loaded in extension only, the stresses for the Current Tie-Bar and the Reduced 
Constraint Tie Bar are identical as expected, with the highest values occurring at the inner edges of 
the outer slots. The maximum stresses for the Elliptical Slot model seem to shift away from the slot 
end slightly. The magnitude is not always lower. This suggests that the elliptical slot design may be 


a benefit to some locations and a detriment to others. Figure IV-5 summarizes these results. 


Extension Only (0.02 inches) B Ellipocat Siot Tie Bar Model 


Current Tie Bar Model 
S Reduced Constraint Tie Bar Mode! 


von Mises Stress (psi) 


267 268 271 272 2-287 288 289 = 290 303 4 305 306 307 308 37 318 
Nastran Element Number 


Figure 1V-5 | Comparison of Critical Element Stresses for Three Designs, Extension 


In the twist only load case, more dramatic results are evident. Again, the Elliptical Slot 
model shifts the stress slightly away from the slot end. It also decreases the maximum von Mises 
stress by nearly 10,000 psi or 17 percent. The Reduced Constraint model shows even more 


spectacular results. This model reduces the maximum von Mises stress by 31,000 psi or 53 
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percent. Not only does the Reduced Constraint model drastically reduce the maximum von Mises 


stress, but it also evens out the stress between the elements. Figure [V-6 summanzes these results. 


‘ . @ Current Tie Bar Model 
Twist Only (0.4 radians) © Elliptcal Siot Tie Bar Model 


§ Reduced Constraint Tie Bar Model 


von Mises Stress (psi) 


267 268 271 272 287 288 289 290 303 304 305 306 307 308 317 318 
Nastran Element Number 


Figure 1V-6 Comparison of Critical Element Stresses for Three Designs, Twist 


The bending only load case causes the models to exhibit behavior similar to that produced 
by the extension load case. Again, the stress for the Elliptical Slot model moves slightly away from 
the slot ends. The stresses produced by the bending load case are only about one-fifth those 


produced by extension and twist. Figure [V-7 summanzes these results. 
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Bending Only (0.006 radians) ih co iin 


S Elliptical Slot Tie Bar Mode! | 
Reduced Constraint Tie Bar Mode! _ | 


von Mises Stress (psi) 
$ 


28 «288280 OB HK S08 07 
Nastran Element Number 


Figure [V-7 Comparison of Critical Element Stresses for Three Designs, Bending 


To summanize, the stresses due to bending are only one-fifth of those due to extension and 
twisting. The bending mode is judged to be not critical, Extension induced stresses at all critical 
locations are within a narrow range. The elliptical slot mainly shifts the maximum stress locations; 
its benefit is not definitely obvious. The twisting induced stresses are dramatically lowered by 50 
percent by the reduced constraint configuration. Serious consideration of this modification is 


warranted. 


Results of the Maneuver 4 flight load case at y = 270 are shown in Figure [V-8. These 
results are formed from a scaled linear combination of the above uncoupled load cases. This 
combined load case indicates that the Elliptical Slot design reduces the maximum stress by about 
5,000 psi or 5 percent. The Reduced Constraint design reduces maximum stress by about 25,000 
psi or 26 percent. Again it is apparent that the Reduced Constraint design evens out the stresses 
among the elements. In fact, the Reduced Constraint design shows an increase in stress for some 
elements, but the maximum stress decreases. It is important to remember that these results are for 


only one flight condition at one rotor position. The trend, however, is significant. 
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Current Tie Bar Model 


Maneuver 4, ‘¥ = 270° 2 Eliptical Slot Tie Bar Mode! 


S Reduced Constraint Tie Bar Model 


100,000 —— = ss 


von Mises Stress (psi) 


287 268 271 272 287 288 289 290 303 304 OS 308 307 308 317 318 
Nastran Element Number 


Figure [V-8 Comparison of Critical Element Stresses for Three Designs, Maneuver 4 


In order to evaluate the cumulative effect of the three designs at the different rotor 
orientations, results were extracted from the finite element analysis. Figure IV-9 and Figure IV-10 
show how the von Mises stress varies for one chosen element as the rotor moves through one 
revolution. The element chosen was the maximum stress element for the flight condition. The solid 
lines are the current design stresses, the dotted lines are the elliptical slot design stresses, and the 
dashed lines are the reduced constraint design stresses. The elliptical slot design lowers both the 
mean stress and stress amplitude slightly. The reduced constraint design lowers both the mean 
stress and stress amplitude significantly for the Maneuver 4 condition and lowers the mean stress 


for the Autorotation 5 condition. 
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von Mises Stress (psi) 


—+— Current Tie Bar, Element 287 —- - & - - Modified Slot Tie Bar, Element 303. = — * — Reduced Constraint Tie Bar, Element 287 


Maneuver 4 


100,000 


90,000 


80,000 


70.000 


40,000 


oO 90 180 270 360 
Rotor Position, Y (degrees) 


von Mises Stress (psi) 


Figure [V-9 Slot End Individual Element von Mises Stresses, Maneuver 4 


Autorotation 5 


te) 180 270 360 
Rotor Position, Y (degrees) 
—#— Current Tie Bar, Element 288 ++ ++Moditied Slot Tie Bar, Element 304 — *— Reduced Constraint Tie Bar, Erement 288 


Figure 1V-10 Slot End Individual Element von Mises Stresses, Autorotation 5 
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It is also interesting to note that the moment required to twist the reduced constraint design 
is significantly lower than that required for both the current and elliptical slot designs. The moment 


required to twist the modified slot design is greatest. The NASTRAN R1 Constraint Moments for 


the three designs with a 0.4 radian twist applied are tabulated in Table [V-1. 


Tie Bar Design R1 (Twist) Constraint Moment* 
Current 
Elliptical Slot 44.0 


Reduced Constraint 13.0 | 


*0.4 radian twist of the two-strap NASTRAN model 
Table [V-1 Tie Bar Constraint Moments 
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Vv. LIFE PREDICTION 


To effectively predict Tie Bar life and evaluate design changes, the probability density 
function, PDF, and the cumulative distribution function, CDF, for the tie bar system are needed. 
Determination of the forms of these distributions requires partitioning the Tie Bar system into its 
components whose distribution functions are known. Conceptually, each Tie Bar strap is first 
partitioned into 4 components, Figure I-6. Each of these components is in turn partitioned into 
smaller volumes or finite elements. In a finite-element model with all elements of comparable size, 
any element could describe a metric volume. The size of this metric is arbitrary; a larger metric can 
be thought of as a combination of smaller metric volumes. A larger metric, however, has a higher 
probability of failure and a lower strength. The selection of the metric volume is such that the 
random vaniable is spatially homogeneous; that is, the stresses are uniformly distributed. The 
uniform stress criterion is the same as selection for element size in finite element analysis; as a 


result, the element can be used directly as a metric. 


A. FLAW DISTRIBUTION 

The mechanistic bases for failure in metals are by flaw growth and by dislocations. For 
brittle failure modes, a given element contains flaws that are either smaller than some critical flaw 
size, a < a, or larger, a 2 a,. Each element can be considered flawed or not flawed based on this 
criterion. In other words, scanning each element for flaw size would create two categories of 
elements, those with flaw sizes greater than or equal to critical, and those with flaw sizes less than 
critical. There are only two possibilities for each element, pass or fail. Respective to the critical 


flaw size, the exceedences are binomually distributed. The binomial distribution has the form: 


f(n) =Cyp"(1- py, 


Equation V-1 


where f(n) is the probability density function, pdf, 
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vy ___! 
” (N=n)Inl’ 


Equation V-2 
p ts the probability that the element has a = a., n is the number of elements that have a > a,, and 


N is the total number of metric elements being considered. In the current investigation, it is the total 


number of elements in each strap. 


When the number of elements is large, and the flaw density is low, then the probability 


density function, pdf, takes on the Poisson distribution form: 


n 


f= "em, 


Equation V-3 
where =Np and is assumed to remain constant. 


For fatigue, the location parameter of the distribution of flaws is tume dependent. As time 
increases, flaws grow and the probability that flaws exceed the critical size increases. The binomial 


distribution then, is a snapshot in time. 


B. COMPONENT LIFE DISTRIBUTION 


For any given instant in time, in order for a given component to have life t, each of its 
elements must have life t. In other words, element | and element 2 and element 3 and ... element n 
must have life t in order that the entire system has life t. Conversely, if any one element does not 
have life t, the whole component fails. In other words, when the weakest-link fails, the system fails. 
In terms of the previous discussion, when any element experiences crack growth to a size greater 


than the critical dimension, the whole Tie Bar component fails. 


In Boolean algebra, the weakest link system reliability is the intersection (74). The 
reliability of the components is conceptually a chain, or series, of links with no parallel or 
redundant links. For a three component system where R denotes reliability or probability of non- 


failure, 
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RA 18%) = RORY, 


X,\PUX,} 


Equation V-4 


Assuming now that the component is made up of many elements and the reliability of each is 


independent, 
R, = REX, VX, 1X40--1.X,} = REX REX REX} REX}, 


Equation V-5 


where R, denotes the Tie Bar component reliability. In this case, the component being discussed is 
itself a system of elements or links in the chain. Failure is assumed homogeneous to mechanism. In 
other words, failure is not related to the size of a component. A larger component will have more 
elements, but the flaw distribution does not change. When this is true, a hazard function =(t) is 


defined which takes on the Weibull form, 
E(r) =r", 


Equation V-6 


where =(r) is the system reliability or hazard function, t is the life of a single element, and a is the 
number of elements in the component. The cumulative distribution function, CDF, can then be 


written as: 
F(t) =1-exp[-=(z)]. 


Equation V-7 


It is important to recognize that the Weibull life distribution contains an underlying Poisson flaw 


distribution, which was described in the previous sub-section. 


C. LIFE CONVOLUTION 


The intrinsic normalized life, t, for a given stress history, S(t), is obtained by convoluting the effect 


of stress via the breakdown rule x. 


39 


t= =f’ K(S(t))at , 


Equation V-8 


where S(t) is a stress norm that is piecewise continuous in time t, fisa non-dimensionalizing and 
normalizing parameter for time, t; is the initial time, tr is the final time, and x is a damage function. 
In the simplest case, this process reduces to the familiar Miner’s Rule, which is a linear 
superposition of damage. 

Different physical processes give rise to different forms of the damage function «(S(¢)). 


The power form and exponential form are explored. Combinations of the two forms are also 


possible. 


1. Power Law Damage Function 


The first proposed damage function is defined using a power relation. This form has been 
observed to fit high cycle fatigue data in metals and is associated with flaw growth. 


Ky b 
xcs =(S2), 


Equation V-9 
where b is a constant exponent, and C, is a constant non-dimensionalizing parameter for stress. 
Both b and C, are material constants. The constants are determined by fitting a line to Strength- 
Life data for a given material in log space where 
1 
~ slope’ 


Equation V-10 


and 
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C, = exp(intercept). 


Equation V-11 


Substituting Equation V-9 into Equation V-8 yields, 


Leu S(t)’ 
raz] (So) at, 


Equation V-12 


and substituting Equation V-12 into Equation V-6 yields, 


Equation V-13 


Finally, substituting Equation V-13 into Equation V-7 yields: 
Ly (Suyy_)” 
fd 
F(t\S) = 1—exp| -| = (50 at| |. 
(«\S) P ; J C 


Equation V-14 


Next, the parameters from the life test are related to the parameters in standard Weibull 


form. Standard Weibull form is: 


R(t) =1- F(t) = exp -(4| 


's 


Equation V-15 


So, Equation V-14 is rearranged and equated to Equation V-15. 
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Equation V-16 


Every element in Equation V-6 is now known. fis an arbitrary normalizing parameter, 
C), b, tg, and a, are determined from material testing, and a is the size effect parameter, which is 
unity if the test specimens are the same size as the actual part. Given a stress history, Equation V- 


6 can be solved for life t. 


2. Exponential Law Damage Function 


The second proposed damage function is defined using an exponential law. This form has 


been observed to fit low cycle fatigue data in metals and is associated with yielding. 


1 S 
K(S(t)) = Jew $) ; 
2 3 


Equation V-17 
where C; is a constant, and C; is a constant non-dimensionalizing parameter for stress. 


Both C, and C; are material constants. The constants are determined by fitting a line to Strength- 


Life data for a given material in semi-log space where: 
C, = slope, 
Equation V-18 


and 


C2 = exp(intercept / C3), 


Equation V-19 


Substituting Equation V-17 into Equation V-8 yields: 


Equation V-20 
and substituting Equation V-20 into Equation V-6 yields: 
\ 4 
= ley lf Sto) 
=(r)=|=| exp) =“ lar) . 
2 a ag 
Equation V-21 


Substituting Equation V-21 into Equation V-7 yields: 


| ( ‘ 
F(t)= iol dex 2) 


3 
Equation V-22 
Next, the parameters from the life test are related to the parameters in standard Weibull form. 


So, Equation V-22 is equated to Equation V-15. 


Equation V-23 


and the equation is solved for t. 


3. Constant Stress Case Example 


Having developed the above, the solution approach for the power law form of the damage 
function is to determine the parameters a, b, and C, from a known stress history. Once known, 


these parameters can be used to predict F(t) of another history. 


For example, for a constant stress case defined as, 


5 0 £5 
tas §, tt," 


Equation V-24 
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S(t) 


Si 


Figure V-1 

b b 
ley S.) =i(4) y 
oH (S) asi(S) Las, 


or 


Equation V-26 


=f 


i 


t 
Equation V-26 shows that life is proportional to elapsed time [ f Fi for a constant 


load S,. Since b is a negative number, an increase in S, reduces life t. 


Now, substituting into Equation V-13, 


(S, ba t ~t, a 
sola) (GY. 


Equation V-27 


and substituting Equation V-27 into Equation V-7, 
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Equation V-28 


which is the cumulative distribution function or cumulative probability of failure for the constant 


stress history. 


Next, the parameters from the life test are related to the parameters in standard Weibull 


form, Equating Equation V-28 with Equation V-15, 


wl) (4) emt -(2) 


Equation V-29 


and the equation is solved for t. 


4. Varying Stress Case Based on Rotor Kinematics 


The Tie Bar stress history is taken to be of the form: 
SVD) = S,, -S, cos y(t) - 90"), 


Equation V-30 


where 


Equation V-32 
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and y(t) is the rotor position in degrees. Figure V-2 shows a typical stress history for the Tie Bar. 


Substituting this new stress history into Equation V-12, 


r= al [S,, - S, cos(w($)- 90°) ae. 


Equation V-33 


Unlike the previous constant load case, the integrand in Equation V-33 is not symbolically 


integrable. Therefore, the integral must either be solved numerically or approximated. 


Tie Bar Stress History 
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Figure V-2 Typical Element Stress History 
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“vi. COUPON TESTS 


Coupon tests were performed by the Tie Bar manufacturer to evaluate different 
manufacturing processes. The test coupons were single, modified production Tie-Bar straps. Each 


coupon had the center two components removed by machining. 


/ “y 
Ye oO 


{O =) 
Ce] | Po) 


e) 


Figure VI-1 — Tie-Bar Coupon 


The following manufacturing process modifications were evaluated either individually or 


in combination: 

e Stamping from new die 

e Pre-twisting to 1800 micro-strain 
e Production glass bead peening 

e High intensity peening 

e Peening at a 45 degree angle 

e Chemical milling 


e Deburming and steel shot peening 


A. STAMPING FROM NEW DIE 


A new die was used to stamp some of the coupons to get a cleaner edge. The original die 


creates a sheared edge that is fairly rough and is thought to provide numerous crack-initiation sites. 


B. PRE-TWISTING TO 1800 MICRO-STRAIN 


Control system checkouts are performed on the ground with the rotors stationary. These 


control checks twist the Tie-Bar approximately 22.5°. Some bending also occurs due to the rotor 
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blade weight. A strain-gaged Tie-Bar installed on an aircraft showed a maximum compression 


strain of -1800 micro-inches during this control check. 


To emulate this check, some coupons (two at a time) were put in a twist test fixture with a 
spacer between them so they would experience the same loading as the outermost straps of an 
entire Tie-Bar assembly. The coupons were twisted to -1800 micro-inches as read from a strain 
gage mounted on the outer strap near a slot end . The coupons were then removed from the twist 


fixture and tested in simple tension-tension until failure. 


C. 6° TWIST LOAD 


Some specimens were subjected to an added load calculated by twisting the coupon 6°. 
Each coupon was not actually twisted. Instead, a strain gage was mounted on a coupon at the slot 
radius and the coupon was mounted in the twist test jig. The coupon was twisted 6° and the strain 
was measured. This strain was converted to stress. The additional coupon load was calculated 
assuming this stress equated to a nominal (P/A) stress. The coupon was then subjected to this twist 


load plus the standard fatigue load. 


The coupons were tested in a tension-tension fatigue machine with a sinusoidal loading. 


The load blocks were repeated until 500 blocks were complete or the coupon failed, except as 
noted. The coupon load blocks are listed in the appendix Table VIII-1. A table of the various 


configurations and their associated lives are also located in the appendix Table VIII-2. 


D. DETERMINING MATERIAL PARAMETERS 


The coupon geometry was modeled using two finite-element models: one for the standard 
coupons, and one for the chemically milled coupons. The chemically milled coupons had 0.002 
inches of material removed from each surface, and the second finite element model reflects this. 
The basic load case for the models was a 0.1 inch displacement on each end of the coupon to put 


the coupon in tension. 
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Figure VI-2 Coupon Finite Element Model 


The T1 Constraint Force was calculated for each model and was scaled for each coupon 
test Load Profile. These scale factors were applied to the finite element models to recreate the 
coupon load conditions. Maximum coupon von Mises stresses were recorded for each Load Profile. 
These maximum von Mises stresses were converted to alternating stresses, S, using the Goodman 


relation. 


Equation VI-1 


where §S, is the fully reversed stress level corresponding to the same life as that obtained with the 
stress conditions, o, and o,,, and o, and o,, are the von Mises stress amplitude and mean stress 
determined from the finite element model. The results are tabulated in the appendix and are plotted 
in Figures V-3 to V-10. These figures have been fitted with lines using a least squares method. In 
the least squares method, each data point is equally weighted and is an expedient method of 


determining the central tendency of the fatigue behavior, i.e. the mean of the S-N curve. 
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Figure VI-6 
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Log-Log S-n Plot of Coupon Test Data, New Stamping, Hi-Intensity Peen 
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Figure VI-7 —__ Log-Log S-n Plot of Coupon Test Data, New Stamping, Chem- Milled 
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Figure VI-8 —Log-Log S-n Plot of Coupon Test Data, New Stamping, Production Peen, 
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Figure VI-9 — Log-Log S-n Plot of Coupon Test Data, New Stamping, Deburred 
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Figure VI-10 Log-Log S-n Plot of Coupon Test Data, New Stamping, Chem-Milled, 
Elliptical Slot 
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For each alternating stress level in Figures V-3 to V-10, there are several values of life. 
There is a mean distribution for each of these stress levels. The fitted lines can be interpreted as the 
mean life of the population at each stress level and may be corrected to the location parameter by 


the relation 
( 
a 


Equation VI-2 


For this conversion, the shape parameter « must be known. To estimate a, a maximum likelihood 
estimator technique can be used to find the parameters of a Weibull distribution that fits the data at 
each stress level. Such extensive data interpretation has not been pursued in the investigation. 
Instead, it was assumed that the strength-life relation is homologous; that is, the fractional life S-N 
curves do not criss-cross. Under this assumption, the shape parameter for strength a, can be 


related to the shape parameter for life a, by 


a 


t > 


wt 
= a, 
Equation VI-3 


where b has already been determined by the least squares fit. Typically, strength data are more 
plentiful and the coefficient of variability cv, in strength are reported in open literature to be 0.5%. 


The coefficient of variation can be converted to the shape parameter by the approximation 


2 


o 


Equation VI-4 


Table VI-1 tabulates the material parameters for the power law development in the 


previous section. Larger values of b indicate longer life. 
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Condition : b C, (psi) 


Production Part, Production Peen 364,171 


New Stamping 2,121,151 


Production Part, 6° Twist 237,081 


143,599 


New Stamping, Hi-intensity Peen 


New Stamping, Chem-milled 168,054 


New Stamping, Production Peen, Chem-milled 195,128 


New Stamping, Deburred -14.03 123,731 
New Stamping, Chem-miulled, Elliptical -27.36 172,509 


Table VI-1 Coupon Test Calculated Parameters 
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vu. CONCLUSIONS 


This research developed a probability modeling method that can aid the prediction of a Tie 
Bar design to meet desired life and reliability levels. The model allows the effects of changes in 
geometry, material, and manufacturing processes to be quantified. Since the model can predict life, 


it could be used for the current Tie Bar to help determine a safe time interval between inspections. 


Two design concepts were explored in this investigation. The elliptical slot design did not 
reduce stress concentration in all flight kinematic configurations. Whether the elliptical slot design 
will yield longer life, is not definitive. The reduced constraint design proposed herein provided 
much more stress concentration reduction as suggested by the finite element study. This reduced 


constraint design may produce large gains in life over the current design. 
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vil. RECOMMENDATIONS 


Some field data has been recorded on current-design Tie Bars. It is recommended that 
detailed breakage data including part hours, number of broken components, and location of broken 
components continue to be recorded. This data can be used to calibrate the model created in this 


investigation. 


While the modeling method developed is basically complete and parts of it have been used 
to calculate results, more effort needs to be spent on organizing the NASTRAN output and feeding 
it through the convolution and probability models. The coupon data could be analyzed more 
thoroughly to provide better material parameters. This investigation only focused on one strap (4 
components) of a two-strap model. Further work should be done to predict the probability of 


failure of at least the eight components of the two outer straps. 


Boeing Vertol is qualifying new designs for the Tie Bar through testing. Results of this 
testing need to be compared to the model results to ensure the model's validity. Field data can be 
compared to both the model results and the testing results to determine what further load and life 


modeling needs to be done. 


Better automation and generalization of the methods descnbed in this investigation would 
be quite useful. A finite element post-processor could theoretically be developed to study the Tie 


Bar as well as any general structure. 
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APPENDIX 


Mean __Load Minimum |Maximum |Block Alternating 
Profile [Load Amplitude |Load Load Cycles Stress 
Ibs Ibs Ibs Ibs psi 
A 1406 940 466 2346 10,000 84,438 
2219 10 87,411 
B | 1406 752 654 2158 10,000 67,550 
1109.5 1109.5 0 2219 10 87,411 
C 1406 470 936 1876 10,000 42,219 
1109.5 1109.5 0 2219 10 87,411 
1876 10,000 42,219 
329] 2304] 2962 132* 70,382 
10,000 35,576 
[| __ 2219] 10 87,411 
1,000,000 44,338 
| 654] 2158 10,000 67,550 
10 
| _ 805] 2007 10,000 
1109.5] ig 2219 10 87,411 
| 1406] —s601/_— 805] ~—S—2007/ 10,000 53,987 
[1481] 1481f 2962 1o| 137,926 
| 466] ~——s*104] 1036] 1,000,000 30,002 
1220] 847|_ Ss 373| S207 10,000 69,935 
961.5]  961.5/ _—iO 1923 10 71,371 
1220[ 847/373] ~~ 2067| ~~ 10,000] ——-6 9,935 


E 

P| 961.5] 961.5] 1923] 10] 71,3711 
IM {| 1220] 1010] 210] 2230] 10,000] 83,3933 
P| 964.5] 961.5] T1923] 0] 71,371] 
IN | 1220] 1174] 46] 2394] 10,000] 96,934 
|| 96.5] 64.5] SS t923] 10] 71,3711 
jo. | 1220] 800] 420] 2020| 10,000] _—66, 0544 


[e) 


a 


1109.5] _ 1109.5] 
* Applied One Time Every Five Blocks 

' 1,000,000 cycles 

? 10,000,000 cycles 

Alternating Stress is calculated using a Goodman relationship and von Mises stress from 


the coupon finite element model. 


Table VIII-1 Coupon Load Conditions 
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Specimen | Production New 6° | Production Deburred| Elliptical | Load Life 
ID Part Stamping | Twist Peen milled Profile (cycles) 
1 | x [| | A 
1Bl_ XxX Et 2 ae eae a 
2B| Xx 2 ees P 
a el ee ee es Re A ; 
Sa ae a ea Stamm 
4 x x =~] B 19,800 
5 xX x | Cc 313,000 
| 2 = Cc 305,000 
7|_ xX 2 D 950,000 
8|_ xX = 2 eee ee D 852,000 
a a ee eS E 1,000,000 
a E 1,000,000 
1i{ xX x 
fax oT CC CT 28.800) 
i} x TT 39.286 
taf xX TT 88,000] 
15] X Pp x fH 28,200] 
ee ed ee 
| | ND a ES aed ees, es ee | aes 
ET) ee 
(| ee 
Se a See ES A (|| 
a ee a eee er ee H 39,700) 
a ee a | es 
| 23st XK 259,700) 
Lat x KT 281.000] 
eS Se” B $8,100 
a a a ee ee ee 183,000 
a OH ht tS 
i i ed H BB08 
a aa a ae a = es 
p80], XT x KT 3,800,000] 
a) ae ee L ci 
a a ee ee ee en 
2! a ae ee eee ee | x | | IH | 4,800,000 
ES) a ae ee 
Lee ik | | ee) 
a a ea ee ee ee ee 
a a ea ee ee 
| a a a ae ee 
i ae i a ee ee ee ee a ee 
ee ee eee eee ee eee aa0G 
a re ee ee | | XT ) pads 
Patt xe oT 3,870,000) 
= eae ee ee eee eee ee M a1.508 
aE ae ee ee ee ee ee L $:000,000 
EC) =: 


an 
N 


Specimen | Production 


New e& 
Stamping | Twist 


Production 


Peen 


High- 45° 


intensity | Peen 


Deburred | Elliptical 


Load 
Profile 


Life 


(cycles) 


R 
Cc 


Table VIII-2 Coupon Test Summary 
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Cc 
M 
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